TAYLOR APPROXIMATIONS OF OPERATOR FUNCTIONS 
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Abstract. This survey on approximations of perturbed operator functions ad- 
dresses recent advances and some of the successful methods. 



1. Introduction 

An active mathematical investigation of perturbed operator functions started in as 
early as 1950's, following a series of physics papers by I. M. Lifshits on the change of 
the free energy of a crystal due to appearance of a small defect. The latter research 
in physics gave birth to the Lifshits-Krein spectral shift function I^TJ |2B] , which 
has become a fundamental object in perturbation problems of mathematical physics. 
Subsequent attempts to include more general perturbations than those in [271 l2H] have 
resulted in consideration of higher order Taylor approximations of perturbed operator 
functions and introduction of Koplienko's higher order spectral shift functions [26| 

Approximation of operator functions also arises in problems of noncommutative 
geometry involving spectral flow (see, e.g., [8]) and spectral action functional (see, 
e.g., [IS]). This investigation was initially carried out independently of the study of 
the spectral shift functions. However, a recent unified approach to the Lifshits-Krein 
spectral function and the spectral flow allowed to establish that these two objects 
essentially coincide [5]. Higher order Taylor formulas have been derived for spectral 
actions in [56], with restrictions on the operators relaxed in [51] by applying more 
universal perturbation theory techniques. 

The proof of existence of the first order (Lifshits-Krein) spectral shift function, 
which is due to M. G. Krein, relied on the theory of analytic functions and was of a 
different nature than the proofs of the other mentioned results on the approximations 
of operator functions. An important object in higher order Taylor approximations is 
the Gateaux derivative of an operator function. When the initial operator and the 
perturbation do not commute, the Gateaux derivative is a complex object, whose 
complexity increases with the order of differentiation. Treatment of such deriva- 
tives and subsequent derivation of Taylor approximations was based on a delicate 
noncommutative analysis, which had been developing for some 60 years. 

To proceed to a detailed discussion of the aforementioned and further results and 
methods, we need to fix some notation. We work with a pair of operators defined in 
a separable Hilbert space H, denoting the initial operator Hq and its perturbation 
V. The perturbation is always a bounded operator and, moreover, some summability 
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restrictions are imposed either on V" or if o- In some instances, Hq is allowed to be 
unbounded, and we will consider only closed densely defined unbounded operators. 
For sufficiently nice scalar functions /, we consider the operator functions f{Ho) and 
/{Hq + V) given by the functional calculus. We are interested in some scalar charac- 
teristics associated with perturbations that are calculated using traces (a canonical 
trace Tr, a Dixmier trace Tr^, a normal trace on a semi-finite von Neumann algebra 
r, and, more generally, any trace rj on a normed ideal X continuous in the ideal 
norm) . 

We consider the remainders of the Taylor approximations 



n— 1 ^ 

Rn,Ho,v{f) = f{Ho + V)-Y,- 



k=0 



f{Ho + tV), 

t=o 



where n G N and the Gateaux derivatives ^|^^g/(-ffo + tV) are evaluated in the 
uniform operator topology. If the n-th order Gateaux derivative is continuous on 
[0, 1], then we have the integral representation for the remainder 

Rn,Ho,v{f) = ^^^3^ (1 - tr-'^f{Ho + tV) dt, (1.1) 

which can be proved by applying functional in the dual space (BiTi))* of the algebra 
of bounded linear operators on "H and reducing the problem to the scalar case. The 
questions we are interested in consist in establishing more specific properties of the 
remainders Rn,Ho,v{f)- 

2. SCHATTEN CLASS PERTURBATIONS 

In this section, we discuss Taylor approximations in the classical setting of pertur- 
bations belonging to the Schatten-von Neumann ideals of compact operators 

S" = {Ae Bin) : := (Tr(|A|"))" < oo}, ae [l,oo) 

(see, e.g., [SU])- The operator functions under consideration come from either poly- 
nomials V or the functions with nice Fourier transforms 

W„ = {/: /(^^/W eLi(R),j=0,...,n}. 

The class W„ includes such widely used sets of functions as C^~^^ (M) and the rational 
functions in Co(M), which we denote by TZ. 

2.1. Spectral shift functions. As a joint finding of many investigations, we have 
the following representation for the Taylor remainders corresponding to self-adjoint 
perturbations of self-adjoint operators. 

Theorem 2.1. If Hq = Hq and V = V* G S*", n E N, then there exists a unique 
real-valued function rjn = rjn,Hoy ^ i^^(M) and a constant c„ > such that 

and 

fiHo + V)-Y,-—^ ^ nHo + tV)\ = Jj^-\t)vn{t)dt, (2.1) 
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for f e Wn. 

The cases n = 1, n = 2, and n > 3 are due to [27] (see also [29]), [26], and 
|42j . respectively. The formula (12. ip has been extended from Wn to the Besov class 
S^j^(M) in [37], [38], and [3], respectively. Differentiability of operator functions in 
the setting most applicable to Theorem 12.11 is discussed in [39] and [6] . The results 
of [HI |39] can also be used to justify that the trace on the left hand side of (12. ip is 
well defined. 

The function ?7„ provided by Theorem 12.11 is called the n-th order spectral shift 
function associated with the pair of operators {Hq,Hq + V). The name to r/i was 
given by M. G. Krein and can be understood from I. M. Lifshits' formula 

ViiX) = TT{EH,Xi-oo, A))) - Tr(EHo+y((-oo, A))) 

holding for Hq and V finite matrices, where Eh denotes the spectral measure of H. 
A number of remarkable connections of the first-order spectral shift function to other 
objects of mathematical physics can be found in the brief survey [9]. More detailed 
discussion of the first order spectral shift function can be found in [121 EDI ET] and of 
the second order one in [23]. When a perturbation V is in the Hilbert-Schmidt class 
S*^, the higher order spectral shift functions can be expressed via the lower order ones 
[T8| [52] . The former are more sensitive to the displacement of the spectrum under 
perturbation, as demonstrated in [SUES]- 
The question of validity of 

Tr(/(i/o + V)- f{Ho)) = [ fit) ri^it) dt, (2.2) 

was also investigated for non-self-adjoint operators Hq and Hq+V. Here the set C C 
is determined by Hq and V. The trace formula (12. 2 p with 1] = T (the unit circle) was 
proved in [2S] for unitary operators Hq and Hq + V such that V E S^. The case of 
arbitrary bounded operators Hq and Hq + V differing hj V E is naturally harder 
than the case of self-adjoint operators. If Hq and Hq + V are contractions, then (12. 2p 
holds with f2 = T and the function rji replaced by a finite complex-valued measure 
I'l, for every / analytic on a disc centered at zero of radius r > 1 [19]. Attempts 
to get more information about the structure of the measure Ui (for example, extract 
an absolutely continuous component) resulted in consideration of only selected pairs 
of contractions and brought to modification of (12. 2p with passage to a more general 
type of integration. The relevant discussion (also for dissipative operators Hq and 
Hq + V) can be found in[Il[21l3ni[Mll3a[3ailZlllHlllS]. 

The higher order version of (12. 2p for pairs of bounded operators has a more plausible 
formulation. 

Theorem 2.2. Let Hq and Hq + V be contractions and assume that V G S"", n > 2. 
Then, there exists a function r]n = Vn,Ho,v L^(T) such that 

f{HQ + y) - ^ _ _ J(Hq + tV)\ = f^-\z) Vn{z) dz, (2.3) 
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for f E V. Furthermore, there exists a constant Cr, 
satisfying (12.31) can he chosen so that 

ii^niii<c„i|y|i::. 



> such that a function ?7„ 



(2.4) 



The case n = 2 for Hq and Hq + V unitaries, where the derivative is evaluated along 
a multipicative path of unitaries instead oi t ^ Hq + tV , is due to |36] (with later 
extension of the class of functions / in [38] ) and for arbitrary pairs of contractions Hq 
and Hq+V joined by the path 1 1— )■ HQ+tV is due to [16] . The case > 3 is established 
in [44] . The spectral shift function r]n satisfying Theorem 12.21 is determined uniquely 
only up to an analytic term (that is, the equivalence class of in the quotient space 
L?-{T)/H^{T) is uniquely determined). Theorem 12.21 can be extended to more general 
functions /. In particular, (12. 3 p with = 2 is established for analytic functions / in 
[T9] , as discussed in Subsection 13.31 for more general traces. 



2.2. Proof strategy. The proofs of Theorems 12. II and 12.21 are very subtle and tech- 
nically involved, so we will give only a flavor of some basic ideas. For simplicity we 
assume that \\Hq\\ < 1, ||iJo + ^11 < V E S'", and f E V. Then our goal is the 
formula 



f^^\z)Ur,{z)dz, 



TT{R^,Ho,v{f)) 

where z/„ is a finite measure, with total variation bounded by 

\\u II < c iii^ir 

II "^n II _ "^n II II n' 

From the integral representation for the remainder (II. ip . we derive 
Thus, if we prove 



(2.5) 



(2.6) 



f{HQ + sV) dt 



sup 

te[o,i] 



n\ \ ds"' 



f{HQ + sV) 



s=t 



<Cn\\V\\:- f 



HI 



(2.7) 



then application of the Hahn-Banach theorem and the Riesz representation theorem 
for the dual space of C(T) implies existence of a measure z/„ satisfying (12. 5p and (12. 6p . 

For 77, = 1, we have 



Tr 



d 

ds 



f{HQ + sV)) =TT{f'{HQ + tV)V), 

s=t / 



which in case of / a polynomial follows from the straightforward calculation of the 
derivative and some combinatorics. Applying the Holder and von Neumann inequal- 
ities then implies (12. 7p with n = 1 and ci = 1. This reasoning does not allow to 
establish the absolute continuity of ui (which was established in [27]), but it can be 
generalized to apply to the higher order case. If, in addition, we take Hq and V to 
be self-adjoint, then application on the spectral theory allows to derive an explicit 
formula for Ui, as it was done in [TO] . 
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Apart from the case of commuting Hq and G 5^, we do not have the convenient 
equahty Tr (^^^^/(i^o + sV)^ = Ti{f"{Ho + tV)V^). However, since the set func- 
tion Ai X A2 ^ Ti [E Ho+tv {Ai)V EHQ+tv{A2)V) , where Ai,A2 are Borel subsets of 
M, uniquely extends to a measure on 



i"^ with total variation HV^Hl, we have 



fiHo + sV) 



(/OW(Ai, A2) ^T{EH,+tv{d\^)VEH,+tv{d\2)V) 



(see, e.g., [53] Theorem 3.12]), which along with the estimate for the divided difference 
IKf )f''lloo < linioo implies (O with n = 2 and C2 = i. 

When n>3, the set function x ■ ■ ■ x A„ t-^ Ti{EHo+tv{Ai)V . . . EHo+tv{An)V) 
can fail to extend to a measure of finite variation on (see [TSl Section 4]). This 
is one of the reasons suggesting that the case n > 3 requires much more delicate 
(noncommutative) analysis of operator derivatives than the case n < 3. 

Pioneering estimates for norms of n-th order operator derivatives are attributed 
to Yu. L. Daleckii and S. G. Krein [17]. In Hq = H* and V = V* e B(n), 
a scalar function / belongs to C^"(M), and the estimates depend on the size of the 
spectrum of the operator Hq. Development of the Birman-Solomyak double operator 
integration (see, e.g., [11]) and subsequent multiple operator integration (see [39] and 
also [0]) resulted in significant improvement of the estimates for operator derivatives. 
It follows from [39] that for Hq = H* and V = V* e S"", 



sup 

tG[0,l] 



Tr 



fiHo + sV) 



< Cn 



\v\ 



11 ' 



where f & B\ 



ool v 



); however, the norm 



is greater than the norm || /^"^ 



The powerful estimates (12. 7p are established in the following theorems. 

Theorem 2.3. ([Sj) If \\Hq\\ < I, \\Hq + V\\ < I, and n e N, then there exists a 
constant c„ > such that for every / G P the following estimates hold. 

(i) If P>n andV e , then 



sup 

t6[0,l] 

(ii) IfV eS"", then 



ds"" 



f{Ho + sV) 



f 



{n)\ 



Il°°(T)" 



sup 

ie[o,i] 



Tr 



ds^ 



fiHo + sV) 



< Cn\\V\ 



\f 



{n)\ 



Il°°(t)' 



Theorem 2.4. ([12]) If Hq = Hq, V = V* , and n eN, then there exists a constant 
On > such that for every f G W„ the following estimates hold. 

(i) If P>n andV e , then 



sup 

te[o,i] 



ds''' 



fiHo + sV) 



<Cn\\V 



f 
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(it) IfV eS"", then 



sup 

te[o,i] 



Tr 



f{Ho + sV) 



<c„||v^llMI/(")| 



L° 



The proofs of Theorems 12.31 and 12.41 (and also analogous estimates for polylinear 
transformations more general than operator derivatives) include a subtle synthesis 
of advanced techniques from harmonic, functional, complex analysis and noncom- 
mutative spaces as well as development of a novel approach to multiple operator 
integration. The principal two cases here are the ones of self-adjoints and unitaries, 
while the case of contractions reduces to the case of unitaries by applying the Sz.- 
Nagy-Foia§ dilation theory 



2.3. Operator Lipschitz functions. Derivation of the estimates of Theorems 12.31 
and 12.41 was preceded by resolution of Krein's conjecture on whether every Lipschitz 
function on M is operator Lipschitz. Detailed discussion of the problem, including 
references to partial results, can be found in [10|H5]: here we only state the concluding 
result and mention some generalizations. 

Theorem 2.5. (@3j) Let f he a Lipschitz function on M. Then, for every a G (1, oo), 
there is a constant c„ > such that 



\\f{B)-f{A)\\^<c^\\B-A 
for all A = A*, B = B* , defined m U with B - A e S' 



Lipi 



a 



The best constant Ca ~ is obtained in [15] . It is known that not every Lipschitz 
function is operator Lipschitz in 5"^ and in BiTi) (i.e., when a e {1, oo}) [20 | |2T | 122]. 
Operator Lipschitzness of functions of normal operators and of functions of several 
variables is discussed in [U [21] . 

3. Some natural generalizations 

If a perturbation V is not compact and no additional restriction on Hq is imposed, 
then the canonical trace Tr of Rn,Ho,vif) is not defined. Depending on the problem, 
one can consider another trace that is defined on Rn,Ho,v{f) for rather general Hq, 
V, and /, or impose extra restrictions on Hq, /, and/or V to ensure Rn,Ho,v{f) ^ "S*^- 

3.1. Compact resolvents and similar conditions. Perturbations that arise in 
the study of differential operators are multiplications by functions defined on M'^, 
which are not compact operators. In this case, the condition V & gets replaced 
by a restriction on the resolvent of the initial operator Hq. 

If Hq equals the negative Laplacian —A and the operator V act as multiplication 
by a real- valued function in L-^(]R^) fl L°°(R^), then 

{Ho- ziy^ -{Ho + V - ziy^ e S\ zeC\R (3.1) 



(see, e.g., [12]). Due to the invariance principle for the first order spectral shift 
function (see, e.g., [12]), the problem for a pair of self-adjoint operators {Ho,V) 
satisfying (13.11) reduces to the problem for a pair of unitaries with difference in S^, 
and ([21]) with n = 1 holds for / G (M) U 7^, as established in [28]. In this case. 
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rji is an element of (M, j^dt) . Existence of the first order spectral shift function 
under more general resolvent conditions is discussed in [251 ESj • 

If i^o = and \^ is a multiplication by a real- valued function in L^(M'^)nL°°(]R'^), 
then instead of the condition (13.11) . we have 

{I + H^)-^/^V e (3.2) 

(see, e.g., [13]). It is established in [26] that for a pair of self- adjoint operators {Hq, V) 
satisfying (13. 2p . there exists 772 G L^(M, j^dt) such that the trace formula (12. ip with 
n = 2 holds for / G 7^. A modified trace formula is obtained in [l3] for a pair {Hq, V) 
satisfying {I + Hq)~^^'^V G S'^. The proofs are based on multiple operator integration 
techniques developed to partly compensate for the lack of the invariance principle 
under the assumption (13. 2p . 

In perturbation problems of noncommutative geometry, typical assumptions on 
the operators are that the resolvent of Hq is compact and V G B{'H). The following 
result is obtained in [51], relaxing assumptions on Hq and V made in [56] . 

Theorem 3.1. Let Hq = Hq he defined in l-L and have compact resolvent and let 
V = V* E B{'H). Let {^k}'k'=i ^6 a sequence of eigenvalues of Hq counting multiplicity 
and let {il)k\k'=i '^'^ orthonormal basis of the respective eigenvectors. Then, for each 
function f G C"+^(]R), with n G N, 

Ti{fiHQ + V))-TT{f{Ho)) 

n— 1 _ 

= E- E (^^n,^..)---(^^V'^n>+Tr(i?^,o,/,„(y)), 

p=l ^ ii,...,ip 

where 

-[i{RH„UV))=0{\\Vr). 

Moreover, the trace formula (12. ip with / G C^(M) is established in [5] for n = 1 
(this is also discussed in the next subsection) and, under the additional assumption 
(/ -|- ifg)~^/^ G S*^, in [51] for n = 2. The respective spectral shift functions r]i and 
ri2 are locally integrable. 

Taylor asymptotic expansions and spectral distributions have also been considered 
in the study of pseudodifferential operators (see, e.g., [T3]). 

3.2. Operators in a semifinite von Neumann algebra. Let be a semifinite 
von Neumann algebra of bounded linear operators defined on "H and let r be a 
semifinite normal faithful trace r on Ai. (The definitions can be found in, e.g., [32].) 
Note that (i3('H),Tr) is one of examples of (A^,r). Let Hq be either an element of 
Ai or an unbounded closed densely defined self-adjoint operator affiliated with Ai 
(that is, all the spectral projections of Hq are elements of A^). The perturbation V 
is taken to be a bounded element of the noncommutative L^-space associated with 
{M, r), that is, 

V e£''= {AeM: := r(|A|")^ < 00}, n G N. 
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Theorem 3.2. If Hq = H* is affiliated with M andV = V* e n e N, then there 
exists a unique real-valued function rjn = rjn^Hoy ^ ^^-^(M) and a constant c„ > such 
that 

WVnWl < Cn\\V\\l 

f{Ho + tV)] = I f^-\t)r^^{t)dt, (3.3) 



and 



k=0 



t=Q 



for f e W„ 



The case n = 1 was established first for a bounded operator Hq in [13] and then 
for an unbounded operator in [7] . The case n = 2 is due to [HI [52] and n > 3 is 
due to [12]. The strategy of the proof is as described in Subsection 12. 2t this strategy 
can be implemented because noncommutative L^-spaces have much in common with 
Schatten ideals (see, e.g., jH]). 

The first order spectral shift function for a pair of r-Fredholm operators differing 
by a r-compact perturbation is known to coincide with the spectral flow [3, Theorem 
3.18]. It is also established in |3] that f l3.3p with n = 1 holds for Hq having r- 
compact resolvent. (In the case (A^, r) = {B{'H), Tr), a r-compact operator is merely 
a compact operator.) 

Theorem 3.3. ([5j) If Hq = Hq is affiliated with Ai and has a r-compact resolvent 
and ifV = V* eM, then, for f e C^{{a, b)), 

T{fiHo + V))=T{f{Ho))+ [ /'(A)r(i?Ho((a,A])-i?^,,+v((a,A]))rfA. 

Analogs of (13.31) with n = 1 and n = 2 for pairs of arbitrary (non-self-adjoint) 
operators in Ai differing by a perturbation V & C"" are obtained in [19]. As to the 
case n > 3, the results of Theorem 12.41 can be extended to pairs of operators in Ai 
by applying dilation of contractions in Ai to unitary operators in semi-finite von 
Neumann algebras constructed in |19j . 



3.3. General traces. The canonical trace Tr is widely used, but it is not the most 
"typical" trace. The distinctive feature of Tr is that it is normal, i.e, has the prop- 
erty of monotonicity. A continuous trace on a normed ideal of compact operators in 
B{T-L) other than 5*^ has a singular component, which vanishes on finite rank opera- 
tors. Detailed discussion of traces and applications of singular traces to classical and 
noncommutative geometry can be found in [32] . 

Let be a semifinite (von Neumann) factor and X a symmetrically normed ideal 
of Ai with norm || ■ \\x- (The definitions can be found, e.g., in [191 122] •) Let tx be 
a trace on X bounded with respect to the ideal norm || ■ \\x- Examples of (X, tx) 
include (S'^,Tr), (£^,r), where r is the normal faithful semifinite trace on A4, and 
{C^^'°°\ Ttuj), where denotes the dual Macaev ideal and Tr^^ the Dixmier trace 

on it corresponding to a generalized limit u on £°°(N). 

The following results are obtained in |19j . 
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Hypotheses 3.4. Consider a set f2, a closed, densely defined operator Hq affiliated 
to Ai, an operator V E I and a space T of functions that satisfy one of the following 
assertions. 

(i) n = conv(a(i/o) U a{Ho + V)) , Hq = H*eM,V = V* , = ^^(R); 

(ii) n = R, Im(i7o) > 0, ln\{HQ + 1^) > 0, and 

= span {A ^ (2 - A)-'= : A: e N, \m{z) < O} ; 

(Hi) Vt = T, \\Hq\\ < 1, \\Hq + V\\ < 1, and T is the set of all functions that are 
analytic on discs centered at and of radius strictly larger than 1 . 

Theorem 3.5. Let Q, Hq, V and T satisfy Hypotheses \3.4\ Then, there exists a 
(countably additive, complex) measure vi = vi^Hoy on Vt such that 

||z/i|| <min{rx(|Re(y)|)+ri(|Im(y)|), \\V\\x] 

and 



Tx{f{Ho + V)-f{Ho))= [ f{X)MdX), 

Jn 



for all f E T . If Hypotheses 3.4^ are satisfied, then the measure Vi is real and 
unique. 

When X = 5^, the measure vi is absolutely continuous, but when X is the dual 
Macaev ideal (with the Dixmier trace), the measure vi can be of any type [T9l The- 
orem 4.4]. Moreover, we do not have an explicit formula for vi in case of a general 
trace tx- Derivation of an explicit formula for vi in case X = S^, Hq = Hq, and 
V = V* relies on the fact that TT(^EHg{-)V) is a (countably-additive) measure, while 
the set function Tt^(^Eho{')^) can fail to be countably-additive (see [T^ Section 3]). 

As another consequence of singularity of Tr^ (and, more generally, of every trace 
satisfying tj(X^) = {0}), we have the following linearization formula. 

Theorem 3.6. Assume Hypotheses \3.4\ and assume tj(X^) = {0}. Then, 

Tx{f{Ho + V)- f{Ho)) = Tx{f\Ho)V). 

Below we consider perturbations in the ideal X^/^ = G Ai : \A\'^ ^-^} 
and impose an additional natural assumption ||Ai?||x < ||y4||ji/2 ||i?||2i/2, which, in 
particular, holds for the ideals S*^, and C^^'^\ 

Hypotheses 3.7. Consider a set Vt, a closed, densely defined operator Hq affiliated 
with Ai, V E X^l"^ and a setT of functions that satisfy one of the following assertions: 

(i) Q = R, lm{Ho) > 0, lm{Ho + 1^) > 0, and 

7 = span {A ^ (2 - A)-'= : keN, lm{z) < O} ; 

(ii) Q = T, \\Hq\\ < 1, \\Hq + V\\ < 1, and T is the set of all functions that are 
analytic on discs centered at and of radius strictly larger than 1 . 

Theorem 3.8. Let fi, Hq, V and T satisfy Hypotheses \3. 7\ Then, there exists a 
(countably additive, complex) measure 1/2 = V2,Ho.v on VL such that 

\\M<\rT{\V?) 
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and 



Ti f{H, + V)- f{H,) - - 



fiHo + tV) 



t=o 



r{\)u2{d\), 



for every f E J^. 
Theorem 3.9. Suppose Txil'^^"^) 



{0}. Either assume Hypotheses 3.1 or else take 



Ho = H* e M,V = V* e X^l^, and = C\W). Then, for every f eT, 



rx { f{Ho + V)- f{Ho) - - 



f{Ho + tV) 



t=Q 



— Tt 

2 ^ 



dt^ 



f{Ho + tV) 



The major components in the proofs of Theorems 13.51 and 13.81 are analogs of the 
estimates (12. 7p . which hold due to the continuity of tx with respect to || ■ \\x- However, 
presence of a singular component in the trace Tx requires more careful treatment of 
the operator derivatives than in the case of the normal trace Tr. 
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